To appear in The Astrophysical Journal 



Scattering by Interstellar Dust Grains. II. X-Rays 

B.T. Draine 

Princeton University Observatory, Peyton Hall, Princeton, NJ 08544; 
draineOastro . princeton . edu 

ABSTRACT 

Scattering and absorption of X-rays by interstellar dust is calculated for a model con- 
sisting of carbonaceous grains and amorphous silicate grains. The calculations employ 
realistic dielectric functions with structure near X-ray absorption edges, with resulting 
features in absorption, scattering, and extinction. 

Differential scattering cross sections are calculated for energies between 0.3 and 10 
keV. The median scattering angle is given as a function of energy, and simple but accu- 
rate approximations are found for the X-ray scattering properties of the dust mixture, 
as well as for the angular distribution of the scattered X-ray halo for dust with simple 
spatial distributions. Observational estimates of the X-ray scattering optical depth are 
compared to model predictions. Observations of X-ray halos to test interstellar dust 
grain models are best carried out using extragalactic point sources. 

Subject headings: dust, extinction - polarization - scattering - ultraviolet: ISM - X- 
rays: ISM 

1. Introduction 

Interstellar grains can absorb and scatter X-rays. The scattering, typically through small 
angles, results in a "halo" of scattered X-rays within ~ 1° of an X-ray source (Overbeds 1965; 
Hayakawa 1970; Martin 1970). The observed properties of these X-ray halos provide a test of 
interstellar grain models. 

The nature of interstellar grains remains uncertain (see Draine 2003a, and references therein). 
This paper will examine X-ray scattering and absorption for a grain model consisting of two separate 
grain populations - carbonaceous grains and silicate grains. With the grains approximated by 
homogeneous spheres with the size distribution found by Weingartner &: Draine (2001; hereafter 
WD01), this grain model is consistent with the observed interstellar extinction law, the observed 
infrared emission from interstellar dust (Li & Draine 2001, 2002), and the X-ray scattering halo 
observed around Nova Cygni 1992 (Draine & Tan 2003). The carbonaceous grains are assumed to 
be graphitic when the grains are large, extending down to very small sizes with the smallest grains 
being individual polycyclic aromatic hydrocarbon molecules. The scattering is dominated by grains 
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with radii a > 100 A, containing > 10 6 atoms; carbonaceous grains in this size range are modelled 
using the optical properties of graphite. 

The primary objective of this paper is to calculate the X-ray scattering and absorption prop- 
erties of this dust model, and to make these results available for comparison with observations. 
The paper has three relatively independent parts, and readers may wish to proceed directly to the 
material of interest to them. 

The first part of this paper - §2 - is concerned with the dielectric functions of graphite and 
MgFeSi04 from infrared to X-ray energies. The resulting dielectric functions satisfy the Kramers- 
Kronig relations as well as the oscillator strength sum rule. 

The second part of the paper - §§3-4 - examines the X-ray scattering properties of interstellar 
dust, including structure near X-ray absorption edges. The angular distribution of the scattered 
X-rays is discussed: the median scattering angle 9 s ^q is found as a function of energy, and a simple 
analytic approximation to the differential scattering cross section is presented. In §4 we provide 
analytic approximations for the intensity distribution in X-ray scattering halos for simple dust 
density distributions. 

The final part of the paper, §5, discusses observations of X-ray scattering halos, and values of 
total scattering cross section derived from these observations. We conclude that the observational 
situation is unclear at this time, due to the combined effects of instrumental limitations and uncer- 
tainties in the spatial distribution of the dust, but the silicate/carbonaceous grain model appears to 
be consistent with the overall body of observational data. Future observations using extragalactic 
point sources could transcend the uncertainties concerning the dust spatial distribution. 

The principal results are summarized in §6 



2. Dielectric Function 

Our objective here is to obtain continuous complex dielectric functions e(u>) = e\ + ie 2 which 
can be used to calculate scattering and absorption by these grains from the submm to hard X-rays. 
For material at LTE (i.e., no population inversions), the dielectric function must have t2 > 0. 
Causality requires that the dielectric function must satisfy the Kramers-Kronig relations (Landau 
& Lifshitz 1960), in particular 

2 „ f°° e 2 (x)xdx 
ei (co) = l + -P 2 \ 1 2 , 1 

where P indicates that the principal value is to be taken. We will proceed by first specifying 62(00) 
at all frequencies, and then obtaining e\(u) using eq. (1). 



In addition to satisfying the Kramers-Kronig relations, the dielectric function must obey the 
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sum rule (Altarelli et al. 1972) 

r , x , 2vr 2 e 2 ^ „ 

/ uj€2{uj)duj = y, n jZj > (2) 

JO m e . 

where nj is the number density and Zj is the number of electrons in element j. Thus we may define 
the effective number of electrons per molecule 

n e e{u) = j X€ 2(x)dx , (3) 

where Vq is the volume per molecule, and n e fr(oo) = the total number of electrons per molecule. 

Ideally, £2(0;) would be measured in the laboratory for the materials of interest. At energies 
below ~20-30 eV the optical constants can be characterized using transmission and ellipsometric 
studies, but calibrated experimental measurements are usually unavailable at X-ray energies. At 
high energies Hco > 100 eV, 62(0;) for a material can be approximated by summing the atomic 
absorption cross sections of the constituent atoms: 



e 2 {u) « - n i E a iA u ) > ( 4 ) 



where aj s is the atomic absorption cross section contributed by electronic shell s of element j. Eq. 
(4) assumes that |e — 1| <C 1, which is valid at X-ray energies. 

At energies well above the photoionization threshold for shell s, atomic photoionization is 
to high momentum free electron states which will have counterparts in the solid, and we can 
approximate a^ s by the atomic photoionization cross section fitting functions a^ % s ((jj) estimated 
for inner shell electrons by Verner & Yakovlev (1995), and for outer-shell electrons by Verner 
et al. (1996), as implemented in the Fortran routine phfit2.f (Verner 1996). Near threshold, 
however, the photoabsorption cross section (jj S {oj) depends on the band structure of the solid, 
leading to "Near Edge X-ray Absorption Fine Structure" (NEXAFS), which could in principle 
permit identification of interstellar grain materials through observations of X-ray absorption and 
scattering near absorption edges (Martin 1970; Woo 1995). 



2.1. Graphite 

We will assume that interstellar grains are constructed primarily from two distinct substances: 
carbonaceous material and amorphous silicate material. In ultrasmall grains, containing less than 
~ 10 5 atoms, the carbonaceous material has the properties of polycyclic aromatic hydrocarbon 
molecules. Because of their small size, scattering by these grains is negligible, and the IR-optical- 
UV absorption can be calculated using absorption cross sections estimated for PAH molecules or 
ions (Draine & Li 2001; Li & Draine 2001). 
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In the larger carbonaceous grains, the nature of the carbon material is less certain. The carbon 
atoms could be arranged in a graphitic structure (pure sp 2 bonding) or there might be a mixture 
of sp 2 (graphitic) and sp s (diamond-like) bonding, perhaps also with some aliphatic (chainlike) 
hydrocarbon material as well. We will assume that the optical response of the carbonaceous material 
in the grains containing > 10 5 atoms can be approximated by graphite, with density p = 2.2 gem -3 . 

Graphite is anisotropic, with the crystal "c axis" normal to the basal plane. In a cartesian 
coordinate system with z || c, the dielectric tensor is diagonal with eigenvalues (ej_, ej_, en). 

The dielectric functions ej_ and en were estimated by Draine &: Lee (1984, hereafter DL84) and 
Laor & Draine (1993). For E < 22 eV we continue to use ey 2 from DL84, while for E > 22 eV we 
use eq. (4) to estimate e|| )2 . For E < 35 eV we use e±^ estimated by DL84, while for E > 35 eV we 
use eq. (4) to estimate e_i_,2- Our final ej_ 5 2 and en 2 are shown in Figures 1 and 3. 



t 1 1 1 1 1 1 — r 




E (eV) 

Fig. 1. — Im(e) for graphite at E < 40 eV (see text). The ex adopted here is equal to ej_ adopted by DL84 below 
30 eV. Also shown, for comparison, are estimates for e± and e\\ by Djurisic & Li (1999). 

Djurisic & Li (1999) recently reestimated e± and en for hw < 40 eV by fitting a model to the 
available experimental evidence. Their estimate for Im(ej_) is shown in Figure 1. We see that it is 
generally similar to the DL84 dielectric function, although with a pronounced peak at ~ 14 eV in 
place of the broader feature peaking at ~ 12 eV in the DL84 estimate of e^. 

For E || c, however, Djurisic & Li concluded that optical ellipsometric measurements are 
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Fig. 2. — Oscillator strength for absorption near the carbon K edge in gaseous anthracene and pyrene (Gordon et 
al. 2003), and for graphite (Shimada et al. 2000, with normalization discussed in text). 



inconsistent with the dielectric function estimated from electron energy loss spectroscopy (EELS). 
Djurisic &z Li obtained eii(a>) from the optical measurements and then independently from the 
EELS data. These two estimates for en - which differ substantially - are shown in Figure 1. The 
Djurisic & Li estimate based on EELS data is fairly similar to the DL84 estimate for en (which was 
based in part on experimental papers using EELS). Given the inconsistencies among the different 
experimental investigations, en in the optical and UV should be regarded as uncertain. 

For carbon, the only X-ray feature is the carbon K edge. We use the measured K edge X-ray 
absorption profile for graphite from Shimada et al. (2000) for 282 - 310 eV. Shimada et al. did not 
determine the absolute absorption strength. We fix the amplitude such that the oscillator strength 
between 282 and 320 eV is 0.38, and we smoothly connect to join the Verner &; Yakovlev (1995) 
photoionization cross section at 320 eV. In Figure 2 we show the oscillator strength density 

df m e V 

dE = W^ £2(tj) (5) 
near the carbon K edge for graphite, and, for comparison, for gas phase anthracene and pyrene 
(Gordon et al. 2003). 

In Figure 3 we show en 2 and ej_.2 from the optical to the X-ray region. Figure 3 shows cim 
and ej_ 1 obtained using eq. (1). 
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Fig. 3. — Im(e) and Re(e - 1) for graphite from 1 eV to 2 keV. 



A useful check on the validity of the dielectric function is provided by eq. (3). Figure 4 shows 
that n e g(u>) w 4 just below the onset of K shell absorption, and n e ^{oj — > oo) — > 6, as expected for 
carbon. 



2.2. Silicate 

We will assume that the silicate material has an olivine composition, Mg2 X Fe2(i_ a; )Si04. Mg 
and Fe are approximately equally abundant in the ISM, and both reside primarily in interstellar 
grains. It is therefore reasonable to take the silicate grain composition to be MgFeSiC^, although 
it is possible that, for example, some of the Fe may be another chemical form. MgFeSi04 olivine 
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has a density 3.8gcm 3 , intermediate between the densities of forsterite (Mg2Si04, 3.27gcm 3 ) 
and fayalite (Fe2Si04, 4.39 gem -3 ), with a molecular volume Vq = 7.5 x 10~ 23 cm 3 . 

For hv < 18 eV we will adopt £2 previously obtained by Draine & Lee (1984) for "astronomical 
silicate" , but with the following modifications to £2 : 

1. the crystalline olivine feature at A -1 = 6.5/im -1 - not seen in interstellar extinction or 
polarization (Kim & Martin 1995) - has been excised with the oscillator strength redistributed 
over frequencies between 8 and 10/um -1 (Weingartner & Draine 2001) 

2. at A > 250/um, £2(0;) has been modified slightly, as described by Li h Draine (2001). For 
250 < A < 1100/xm the revised e 2 is within ±12% of e 2 adopted by DL84. 

Above 30 eV, we use eq. (4) to estimate €2 from atomic photoabsorption cross sections, except 
near absorption edges (see below). Between 18 and 30 eV, £2 is chosen to provide a smooth join 
between £2 from DL84 and £2 estimated from the atomic photoabsorption cross sections. 

The threshold energy for photoionization from the K shell of atomic oxygen is 544.0 eV for 
ionization to O II ls2s 2 2p 4 ( 4 P), and 548.9 eV for ionization to O II ls2s 2 2p 4 ( 2 P); the strong Is — 2p 
absorption line, with FWHM ss 0.14 eV, lies at 527.0 eV (Stolte et al. 1997). From the theoretical 
line profile of McLaughlin & Kirby (1998), the Is — 2p transition has an oscillator strength / ~ 0.10. 

X-ray spectroscopy of several galactic X-ray sources using the Chandra X-Ray Observatory has 
detected a strong and narrow absorption line at 527.5 eV which must be the O I Is — 2p transition, 
and a nearby absorption feature at 530.8 eV, with FWHM ~1.0eV (Paerels et al. 2001; Schulz 
et al. 2002; Takei et al. 2002). Paerels et al. and Schulz et al. suggest that the 530.8 eV feature 
is due to iron oxides, possibly Fe2C>3. However, O II is expected to have Is — 2p absorption at 
approximately this energy, and the interstellar O II abundance is probably large enough that the 
resulting absorption feature should be conspicuous. It therefore seems likely that the observed 
narrow feature at 530.8 eV is (at least primarily) O II Is — 2p. 

In the absence of published X-ray absorption spectra for amorphous silicates, we estimate the 
oxygen K edge absorption using measurements on crystalline olivines. Li et al. (1995) have measured 
Mg K edge and Si K edge absorption for forsterite Mg2SiC>4, and Henderson et al. (1995) have 
measured Fe K edge absorption in fayalite Fe2Si04. We adopt these profiles for the corresponding 
K edges in amorphous olivine MgFeSiC>4, with the absorption profile strengths adjusted to match 
the Verner et al. photoionization cross sections for the corresponding atoms well above threshold. 

Fe L2,3 edge spectra for a number of minerals have been studied recently by van Aken & 
Liebscher (2002). The L2 edge corresponds to the remaining 2p 5 electrons in a 2 P3/2 term, and the 
L3 edge when they are in a 2 Pi/2 term. The Fe in olivine is 100% Fe 2+ . The spin-orbit splitting for 
Fe 2+ produces a 12.8 eV separation between the L2 and L3 maxima (van Aken &, Liebscher 2002), 
with the L2 peak at 721.3 eV and the L3 peak at 708.1 eV. Following van Aken & Liebscher, we 
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model the near-edge absorption in olivine by 



2,3 



00 



3vr 



arctan( — (E — E a ) + — 

w n , / 2 



+ 



3vr 



arctanf — {E - E b ) J + J 



+ ^ exp 



i=l 



£7- £7,- 



(6) 



with £7 a = 708.65 eV, Ef, = 721.65 eV, and w a = Wb = leV. The gaussian components are 
located at E x = 707.8 eV, E 2 = 710.5 eV, E 3 = 720.6 eV, and £ 4 = 723.3 eV. van Aken & 
Liebscher do not give values for the widths of the gaussians, but inspection of their Fig. 1 suggests 
Wl = w 2 = w 3 = w 4 = 1.25 eV. Minerals with 100% Fe 2+ have Ai/A 2 « A3/A4 « 5, and 
A3/A1 ~ 0.28. For 100% ferrous Fe we estimate A± ~ 6 from the spectrum of ilmenite FeTi03 in 
van Aken, Liebscher & Styrsa (1998). To determine the normalization <7o, we require that eq. (6) 
match the Verner et al. cross section at 800 eV. 

For some absorption edges no measurements are available for appropriate minerals. Since the 
near-edge absorption is proportional to the density of electronic states above the Fermi surface, 
the near edge spectra for different elements in a compound show considerable similarity in the 
dependence on energy, but with displacements in energy due to the different binding energies of 
the electron being excited. We use the Si K edge in forsterite Mg2SiC>4 (Li et al. 1995), shifted in 
energy by AE = — 1308 eV [the difference in K edge ionization thresholds for Si I (1846 eV) and 
O I (538 eV)] to determine the O K edge absorption profile up to 575 eV. This results in onset of 
O Is absorption at ~527.8 eV, and an O Is absorption peak at 537.6 eV. The absolute absorption 
coefficient is fixed by requiring that the absorption at 575 eV match that calculated from the atomic 
photoionization cross sections. 1 

A similar approach is taken for the other absorption edges, using measured profiles for Mg K 
edge and Si K edge absorption in forsterite Mg2Si04 (Li et al. 1995), and Fe K edge absorption in 
fayalite Fe2Si04 (Henderson et al. 1995) (see Table 1). 

Although this procedure is not expected to provide accurate estimates of the near-edge ab- 
sorption spectra for amorphous silicates, the provisional near-edge absorption profiles so obtained 
will provide a realistic example of the kind of near-edge absorption and scattering expected from 
interstellar grains. 

Figure 5 shows Re(e) obtained from Im(e) using eq. (1). Figure 4 provides a check on the 
adopted dielectric function: lim^^oo n e fr = 84.01, in agreement with the value of 84 expected 
theoretically for MgFeSiC>4. 



Absorption near the O K edge in Si02 has been measured by Marcelli et al. (1985), who find an absorption peak 
at 540 eV in both a— quartz and glassy Si02- EELS studies by Wu et al. (1996) and Sharp et al. (1996) show the 
peak at 535 eV for quartz, while Garvie et al. (2000) use EELS to locate the quartz peak at 537.5 eV, with the 
absorption edge located at ~ 535 eV. The Si02 polymorphs a-quartz, coesite, and stishovite have their absorption 
peaks within ~ 1 eV of one another (Wu et al. 1996), and likewise their absorption edges agree to within ~ 1 eV. 
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Table 1: X-ray Edge Absorption Parameters 
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"Energy at onset of absorption 
^Energy at peak absorption 

Teak absorption cross section/atom contributed by this shell. 

d Energy shift relative to adopted profile. 

e Shimada et al. (2000) 

/ Li et al. (1995) 

9 Henderson et al. (1995) 

h van Aken & Liebscher (2002) 

3. Cross Sections for X-Ray Absorption and Scattering 

Adopting the dielectric functions discussed in §2, we calculate the scattering and absorption for 
the WD01 dust grain mixture. We use the Mie scattering theory program of Wiscombe (1980, 1996) 
for x = 2-ira/X < 2 x 10 4 , and anomalous diffraction theory (van de Hulst 1957) for x > 2 x 10 4 . 
Anomalous diffraction theory provides an excellent approximation at the X-ray energies where 
x = 5067(a/^m)( J B/keV) > 2 x 10 4 (see Figure 7 of Draine & Tan 2003). 

The scattering and extinction cross section calculated for the dust mixture is shown in Fig. 
6 for 0.1-10 keV X-rays, with the six strongest absorption edges shown in Fig. 7. Also shown is 
the absorption per H atom due to gas-phase absorption calculated using phf it2.f (Verner 1996), 
with interstellar gas-phase abundances. At energies 13.6 < E < 250 eV, absorption by neutral H 
and He is very strong, making it difficult to observe the dust absorption and scattering. Above 250 
eV, however, observations of extinction and scattering by dust become feasible for suitably bright 
sources on sightlines with sufficient dust columns. At E > 800 eV the extinction is primarily due 
to dust grains. 
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Fig. 4. — Effective number of electrons n e g as a function of photon energy for graphite and MgFeSiCU. 

As seen in Figure 7, the calculated scattering cross sections show conspicuous structure in the 
vicinity of the major absorption edges. This occurs because at X-ray energies Re(e — 1) tends to 
be negative, and an absorption feature increases Re(e) (reducing |e — 1|) just below the absorption 
feature, and decreases Re(e) (increasing |e — 1|) just above the feature. Since the scattering is 
approximately proportional to |e — 1| 2 , this results in a reduction in scattering below an absorption 
feature, and an increase in scattering above. 2 This argument applies to 5 of the 6 absorption edges 
in Fig. 7; the exception is the C K edge, for which e — 1 actually becomes positive below the 



2 Takei et al. (2002) estimate that the dust scattering cross section would be reduced at energies just above the 
O K edge. We find, to the contrary, that the dust scattering cross section is increased just above the O K edge - see 
Fig. 7. 
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Fig. 5. — Im(e) and Re(e - 1) for silicate MgFeSi0 4 . 

absorption edge, with a local peak in scattering just below the K edge. 

For the lower-energy absorption edges (C K, O K, Fe 1^,3) there is significant variation in the 
scattering cross section near the absorption edge. As a result, the extinction profile is not the same 
as the absorption profile. In the case of the O K edge, the extinction peak is at 538.0 eV, whereas 
the absorption peak is at 537.6 eV. The energy-dependence of the scattering optical depth can be 
determined by dividing the spectrum of the scattered X-rays by the spectrum of the point source 
component. If this can be done with a signal-to-noise ratio > 10 for ~ 2 eV bins, one should observe 
the structure in <r sca seen near 285, 538, and 708 eV in Figure 7. 

At X-ray energies, the dielectric functions of grain materials become close to unity (see Figures 
3 and 5), the wavelength is small compared to the typical grain radius, and the grains are very 
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E (keV) 

Fig. 6. — X-ray extinction and scattering cross section per H nucleon due to dust, and absorption due to gas. Data 
available at http://www.astro.princeton.edu/~draine/dust.html 



strongly forward-scattering. The quantity 

dcr sca da da 2 da 

di^ = n dn =2 * {1 - cose) dn* 7T9 dn (7) 

is proportional to the number of scattered photons per logarithmic interval of scattering angle; the 
location of the peak shows the "typical" scattering angle. Figure 8 shows ixQ 2 da jd£l for the WD01 
dust mixture for selected energies from 0.3 keV to 10 keV. 

Let 9 s ^o(E) be the median scattering angle for photons of energy E. Figure 9 shows 6 s ,5o(E) 
for the WD01 dust mixture; comparison with Fig. 8 shows that, as expected, the median scattering 
angle # s ,5o is very nearly the same as the angle where ir9 2 da/d£l peaks. Also shown in Fig. 9 are the 
10th and 90th percentile scattering angles, 6 s ,io and s ^q, corresponding to 10% and 90% enclosed 
power. Note that to a very good approximation, 

Qs,io ~ > #8,90 ~ 36^50 • (8) 

For E > 0.5 keV, the median scattering angle for the WD01 dust mixture can be approximated by 

S , 5O « 360" (^pj . (9) 
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Fig. 7. — X-ray extinction and scattering cross section per H nucleon near major absorption edges. 

The median scattering angle for a circular aperture of diameter d is 0.53A/<i (Born & Wolf 1999), so 
equation (9) corresponds to the median scattering angle for an aperture of radius 0.19/iin, consistent 
with the size of the grains which dominate the visual extinction and polarization of starlight, account 
for most of the interstellar grain mass, and are expected to dominate the X-ray scattering. The 
smaller grains, while more numerous, make only a minor contribution to scattering at E > 0.5 keV. 
Note that at lower energies, the median scattering angle rises above the approximation (9), due in 
part to the increasing importance of smaller grains (which contribute most of the geometric cross 
section of the grain population) at these energies. 



For this dust model, the differential scattering cross section can be approximated by the simple 
analytic form 
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Fig. 8. — 2ir8 2 da/dn versus scattering angle 9 S at selected energies. This function peaks at approximately the 
median scattering angle. The broken curves show eq. (11) for E = 1.0 and 5.0 keV. 



with the total cross section for scattering angles < 6 



0"sca(< 9) = cr s 



(9/9 



s,50j 



l + (9/9 



(11) 



Eq. (11) reproduces the empirical result that # Sj io = # s ,5o/3, and 9 s> go = 36> Sj 5o. The approximation 
(10) is plotted in Fig. 8 for E = 1.0 and 5.0 keV, showing that it does indeed provide a good fit. 



4. X-Ray Scattering Halos: Models 
4.1. Models 

For a point source at distance D, scattering by dust on the sightline a distance r = xD from 
the observer produces a scattered halo around the point source (see, e.g., Draine & Tan 2003) with 
the halo angle 9^ related to the scattering angle 9 S through 



9 h « (1 - x) 



(12) 
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Fig. 9. — Median scattering angle # s ,5o as a function of energy for the WD01 grain model. Also shown are scattering 
angles 6 a ,io and # s ,9o for 10% and 90% enclosed power. Broken lines show asymptotic behavior for E > 1 keV. 



Let iVhaio be the total flux of singly-scattered photons, and A r v ia i Q (< Oh) be the flux of photons at 
halo angles < Oh- Define the fraction of halo photons interior to Oh'- 



g(0 h ) 



(13) 



If the dust density is assumed to be plane-parallel perpendicular to the sightline, then for the 
small-angle scattering appropriate to X-ray energies, the scattering halo is given by 



9{Vh) 



where the dimensionless dust density 



p(x) 



p(xD) 



(14) 



(15) 



J p(xD)dx 

where p(r) is the dust density along the sightline at distance r from the observer. 

Because the differential scattering cross section for the WD01 dust mixture can be approxi- 
mated by eq. (10,11), we have 



p(x) [l + (l-x) 2 (0 s , 50 /0 h f 



(16) 
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If the scattering is by a single sheet of dust at distance x^D, then p{x) = S(x — x^) and 



{9 h /e h , 



g{o h ) « 1 r ( ' R r :f' v ' ^,50 =(1-^,50 . (17) 

1 + (Vh/Uh,50) 

For a uniform dust density gradient (with (3 = corresponding to uniformly-distributed dust) 

p(x) = (l-0) + 20x (-1</9<1) , (18) 
eq. (14) can be integrated to obtain 

g(6 h ) = (l + (3)^ a rct a n(0s, 50 /0 h )-p(-^) In [l + (6 8t50 /e h ) 2 ] ; (19) 

g{6h) is plotted in Fig. 10 for 5 cases: (3 = —I (1/4 of the dust between x = 0.5 and 1); [3 = —.5 
(3/8 of the dust between x = 0.5 and 1); (3 = (uniform dust); (3 = 0.5 (5/8 of the dust between 
x = 0.5 and 1); [3 = 1 (3/4 of the dust between x = 0.5 and 1). Also plotted is the case where the 
dust is all at x <C 1, with 9^ = 6 S . For the above dust distributions, Table 2 gives the halo angles 
0fc,iO! 0ft,5O> 0/i,9O enclosing 10%, 50%, and 90% of the halo power for single-scattering. 
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Table 2: Halo Structure Parameters for WD01 Dust 



dust density distribution 
/3 = -1 (3= -0.5 uniform (3 = 0.5 (3 = 1 



@h,w/9s,50 


0.1663 


0.1032 


.0664 


.0467 


.0353 




0.631 


0.530 


0.429 


0.337 


0.262 


#h,9o/#s,50 


2.084 


1.882 


1.660 


1.413 


1.139 


MO x (S/keV) 


59.9" 


37.2" 


23.9" 


16.8" 


12.7" 


6> M0 x (S/keV) 


227" 


191" 


154" 


121" 


94.3" 


M o x [E/ keV) 


750" 


678" 


598" 


509" 


410" 



5. X-Ray Scattering Halos: Observations 

The total cross section for X-ray scattering can be measured by imaging the scattered X-ray 
halo. The flux of scattered photons, iVh a io, is related to the flux in the point-source component, 

Nptsrc, by iVptsrc = (^ptsrc + ^halo) exp(-T sca ), SO 

r sca = ln(l + AWAWc) = - In (1 - /halo) , (20) 

^ ^ -^halo _ (21) 

A^halo + N ptslc 

cq. (20) is valid even when multiple scattering takes place. Absorption (by dust or gas) has a 
negligible effect on Ar ha i /A^ ptsrc , because unscattered and scattered photons are affected essentially 
equally. Estimation of r sca from (20) requires determination of the flux of scattered photons AT ha i 
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Fig. 10. — Fraction of the single-scattering halo falling within halo angle 8^, as a function of 9h/9 a ,m, where 6 a ,s,o 
is the median scattering angle. 



integrated over all halo angles. If the halo flux is measured only for halo angles 0\ < 0h < 02, the 
total halo flux can be estimated from 

JVhalo = T^T , (22) 

g(0 2 ) -gipi) 

the function g(6h) depends, of course, on assumptions concerning both the grain model and the 
distribution of dust along the sightline. 



5.1. Cen X-3 

Woo et al. (1994) used the ASCA X-ray observatory to measure the X-ray halo toward the 
massive X-ray binary Cen X-3 (Ay rs 4.3) as a function of orbital phase, at 1.5 and 2.5 keV. Their 
results for r sc& /Ay are shown in Figure 11. Also shown is the ratio T sc& /Ay estimated for the 
WD01 grain model. Woo et al.'s values of T sca _/Ay at 1.5 and 2.5 keV are a factor ~3.5 below the 
prediction of the WD01 grain model. Can we understand this? 

At 1.5 keV, the ASCA point spread function was such that even when the point source compo- 
nent was at a minimum, the halo intensity exceeded the point source intensity only for 6^ > 200"; 
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Fig. 11. — Scattering optical depth T sca per unit visual extinction Av as calculated for the WD01 model (solid line) 
and as estimated from various observations (see text). Chandra imaging of GX 13+1 at 50" < 6 S < 600" (Smith, 
Edgar, & Shafer 2002) has been corrected for photons interior to 50" assuming uniformly distributed dust (open 
squares) or dust with density proportional to distance (filled squares). Chandra observations of Cyg X-l at 8 S < 120" 
(Yao et al. 2003) are uncorrected for photons beyond 120" (open diamonds) or corrected assuming uniform dust 
(filled diamonds). See text. 



when the point source was at maximum the halo intensity exceeded the point source intensity only 
f° r @h ^ 450". For uniformly-distributed dust and E = 1.5 keV, we expect a median halo angle 
&h,m ~ 100" (see Table 2); it is therefore clear that the ASCA observations were insensitive to 
most of the scattered photons. It therefore seems plausible that the scattered flux may have been 
underestimated by factors of 2-4 due to the dominance of the point source profile at halo angles 
9 < 200 — 400". The determination of r sca /Ay by Woo et al. should be treated as a lower bound 
rather than a measurement. 



5.2. ROSAT Observations 

ROSAT observed both unscattered and scattered X-rays from Nova Cygni 1992 at a number 
of epochs (Krautter et al. 1996). The most recent reanalysis of the ROSAT data found r sca = 
0.211 ± 0.006 at the median photon energy ~480 eV (Draine & Tan 2003). Adopting E(B - V) = 
0.19 estimated from observations of Ha/H/3 (Barger et al. 1993; Mathis et al. 1995), we obtain 
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T sc& /Ay = 0.36. As seen in Figure 11, this value is in excellent agreement with the WD01 grain 
model. Draine & Tan have also carried out detailed modelling, and conclude that the observed 
X-ray halo profiles at 9 different epochs are in good agreement with the WD01 dust model. 

Predehl & Schmitt (1995, hereafter PS95) used ROS AT to estimate / ha i for Cyg X-l, GX 13+1, 
and 19 other galactic sources for which Ay was also available. In each case, PS95 have fitted the 
observations with the ROS AT psf plus a theoretical dust model, and used this to estimate the total 
number of scattered photons. Figure 11 shows the resulting r SCSL /Ay versus the average photon 
energy for each source. The estimates of T sca /Ay inferred from the PS95 observations are generally 
a factor ~ 2 — 4 below the WD01 grain model. If the PS95 values of /halo are accurate, this would 
indicate a serious problem with the WD01 grain model. 

However, for the 21 sources, the halo angle at which the intensities of the fitted halo and psf 
were equal was 70" or larger; for 50% of the sources this angle was 130" or larger. At a typical energy 
of ~1.2 keV, the median scattering angle is 9h 50 ~ 130" for uniformly-distributed dust (see Table 
2). Estimates of /halo therefore rely heavily on the dust model to separate the halo from the point 
source at small halo angles. The modelling by PS95 employed power-law grain size distributions 
dn/da oc a~ q for a < a max - For the 21 sources, the median q value was 4.0 and the median value 
of o max was 0.18/xm. Power-law grain size distributions with q « 3.5 and a max ~ 0.25/xm (Mathis, 
Rumpl, & Nordsieck 1977; Draine & Lee 1984) provide a good fit to the interstellar extinction, but 
the size distributions adopted by PS95 - with generally steeper power laws and smaller values of 
Omax _ had insufficient mass in large grains 3 and do not provide a good fit to the extinction law. 
The scattering at small halo angles is dominated by the larger grains, so the PS95 model-fitting 
may have systematically underestimated the actual halo intensity at small halo angles. 

Given likely uncertainties in the psf and the model-fitting, it seems likely that the PS95 values 
of T sca /Av are systematically low, perhaps by factors as large as 2-4. We now examine two particular 
sightlines. 

5.3. GX 13+1 

GX 13+1 is a low-mass X-ray binary system at / = 13.5°, b = 0.10°, and an estimated distance 
D ~ 7kpc, corresponding to a distance above the plane Dsinfe ~ 12 pc. Garcia et al. (1992) 
estimated Ay < 14.4. X-ray spectroscopy gives iVn = 2.9 ± 0.1 x 10 22 cm~ 2 (Ueda et al. 2001), 
corresponding to Ay = 15.5 for the standard conversion Nu/Ay = 1.87 x 10 21 cm~ 2 (Bohlin, 
Savage, & Drake 1978). We adopt Ay = 14. The mean photon energy for the ROSAT observations 
is E= 1.69 keV (PS95). 

PS95 find a halo fraction /h a io = 0.335 for GX 13+1, but (as discussed above) this may be an 



3 The median q = 4.0 corresponds to equal mass per logarithmic interval; for q > 4.0 the grain mass diverges at 
small radii unless a lower cutoff is imposed. 
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underestimate: 



• For GX 13+1, PS95 estimate that the halo and psf have equal intensities at ~75" (see their 
Fig. 10). 4 For uniformly-distributed dust and E = 1.69 keV, the median halo angle dh,50 ~ 
91", but on this sightline the dust density may be enhanced closer to the source, leading to 
a reduction in 9h,5o', f° r a linear gradient with (3 = 1, 6h,5o ~ 58". Thus it appears possible 
that > 50% of the scattered photons may have been misattributed to the psf. 5 

• For GX 13+1, PS95 used a dust model with dn/da oc a~ 3 - 8 and a max = 0.18/xm; as discussed 
above, this underestimates the abundances of a > 0.1/tm grains, and therefore underestimates 
the contribution of scattering at small halo angles. 



It therefore appears that PS95 could have underestimated iVhaio- While difficult to quantify, it 
seems possible that the true value of /halo might be as large as 0.65 (the value predicted by the 
WD01 grain model for E = 1.69 keV and Ay = 14). 

The X-ray halo around GX 13+1 has recently been observed by the Chandra X-ray telescope 
(Smith, Edgar, & Schafer 2002) at energies between 2.1 and 3.9 keV, with energy resolution ~ 
0.2 keV. Phenomena referred to as "pileup" and "grade migration" in the detection system affect 
the Chandra ACIS images as far as 50" from the source. Smith et al. also discuss the current 
uncertainties concerning the Chandra psf at angles > 50". Using a preliminary psf based on 
observations of Her X-l, they estimate what they refer to as "total observed halo fraction" 1(E) 
by integrating the psf-subtracted and background-subtracted count rates from 50" to 600", and 
dividing by the estimated psf count rate in the absence of saturation effects. Thus 

I{E) _ (l-e-^MmO")-g(50")] ^ (23) 

If we assume a model for the dust distribution, we can use g(9h) to obtain 

I 



In 



1 + 



[5(600") -5(50")] J 



(24) 



where for a uniform gradient g(9h) is given by eq. (19). In Figure 11 we show the values of 
T sca /Ay obtained from the observed I, with g(6h) calculated for the WD01 grain model assuming 
(a) uniformly-distributed dust ((3 = 0) and (b) dust with a density gradient (3 = 1. Given the 



4 The ROSAT psf fit given by Boese (2000) has 90%-enclosed-power radii of 32" at 0.5 keV, 25.5" at 1 keV, 35" at 
1.5 keV, 52" at 1.69 keV, and 95" at 2.0 keV. It is not clear how close the actual psf is to the fit given by Boese. 

5 The ROSAT imaging extends to 8 h w 2000", but for GX 13+1 the background is estimated to exceed the halo 
intensity for du > 900"; determination of the background is itself difficult, and underestimation of iVhaio due to 
background oversubtraction at 9 > 600" is an additional possibility. However, since 90% of the scattering at 1.69 keV 
is at scattering angles < 1080"/ 1.69 = 640", underestimation of the halo intensity at Oh > 600" would have only a 
small effect on iV h aio- 
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location of this source, we expect the dust density to be increasing toward the source, so the (3=1 
model is reasonable. For (3 = 1, the inferred r sc& /Ay is a factor ~1.5 below the WD01 model 
for 2.4-3.4 keV, and a factor 1.1-1.3 below the WD01 model at 3.4-4.0 keV. These observations 
suggest that the WD01 model may overestimate T sca /Ay, although it should be kept in mind that 
these results required substantial corrections for unobserved halo interior to 50". 

The best way to use the information in the oberved scattered halo is to try to reproduce the 
observed radial profile of the scattered halo using a dust model, and Smith et al. tested various 
grain models in this way. For dust distributed uniformly between source and observer, for the 
WD01 grain model they find a best-fit gas column Nr = 1-65 x 10 22 cm~ 2 , significantly smaller 
than the value 2.9 ±0.1 x 10 22 cm -2 estimated from ASCA observations (Ueda et al. 2001). Note, of 
course, that additional dust could be located at x > 0.75 without appreciably affecting the observed 
I(9h > 50"), since the additional halo contribution would be mainly below the 50" lower cutoff. 6 



5.4. Cygnus X-l 

Cygnus X-l consists of an O star primary with a black hole companion (see Tanaka & Lewin 
1995), located at I = 71.33°, b = 3.07°, and an estimated distance D = 2.5 ± 0.4kpc (Bregman 
et al. 1973; Ninkov et al. 1986), placing it at a height Dsvub « 130 ± 20 pc above the plane. 
The 09.7 lab primary is reddened by E(B — V) = 1.12 (Bregman et al. 1973), corresponding to 
Ay = 3.5; this is consistent with N-r = 6.2 x 10 21 cm~ 2 from X-ray absorption spectroscopy (Schulz 
et al. 2002). Based on studies of reddening vs. distance for stars within 50' of Cyg X-l (Bregman 
et al. 1973; Margon et al. 1973) it appears that the dust is distributed approximately uniformly 
along the sightline. 

PS95 observed Cygnus X-l with ROSAT, and found /halo = 0.116 at 1.2 keV. The psf and 
scattered halo intensity were estimated to be equal at w 110". The PS95 result r SCSL /Ay = 
.037 mag" 1 at 1.2 keV is a factor 3 below the WD01 model (see Fig. 11). 

Yao et al. have recently used Chandra observations of Cygnus X-l to infer the scattered halo, 
using a technique designed to minimize the effects of "pileup", allowing the excellent angular 
resolution of Chandra to be used to observe at small halo angles. Yao et al. neglected halo angles 
6 > 120", but were able to measure the halo as close as 1" from the point source. Their fractional 
halo intensity FHI is the ratio of the halo counts interior to 120" divided by the counts from the 
psf plus the halo interior to 120", and is related to r sca by 

FHI — - . (25) 

Yao et al. also report the radius 0* containing 50% of the halo counts within 120" of the source, i.e., 
<?(#*) = 0.5(7(120"). In Figure 12 we show the variation of 6* with E calculated for: (1) uniform 



6 0h,m = 43" for E = 2.1 keV and dust at x d = 0.75. 
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dust; (2) thin sheet at Xd = 0.7; (3) thin sheet at Xd = 0.8; (4) thin sheet at Xd = 0.9. The single 
sheet models clearly do not reproduce the observations. The uniform dust model is in approximate 
overall agreement with the distribution of halo counts within 120", although it does not reproduce 
the concentration of the halo at E < 3 keV found by Yao et al. 



i 1 — i — r 




E (keV) 

Fig. 12. — Radius 9* enclosing 50% of the halo counts within 120" of the point source. Data points are observations 
of Cyg X-l by Yao et al. (2003). 

Adopting the uniform dust model, we can now use the WD01 model to correct for halo counts 
at 9 h > 120": r sca = ln{l + FHI/[(1 - FHI) 5 (120")]}, with the results plotted in Fig. 11. At 
E < 2.5 keV, where the estimated fractional uncertainties are smallest, the values of t SCSl /Av found 
from the Yao et al. results fall a factor ~ 1.5 — 2 below the predictions of the WD01 model. The 
reason for this is not apparent. Perhaps the WD01 model has overestimated T SCSu /Ay; alternatively, 
perhaps the background has been overestimated - background oversubtraction would be consistent 
with the surprising concentration of the halo seen at E < 2.5 keV in Fig. 12. 

For 3-4 keV, the inferred values of T sca /Av, and the measured halo half-light radii in Fig. 12, 
are consistent with the WD01 dust model. However, at 5 — 7keV the values of T sca /Ay measured by 
Yao et al. exceed the predictions of the WD01 model, although the error bars are now large. It is 
difficult to envision a dust model which could have such a large value of t SCSl /Av at these energies, 
so we suspect that the observations are affected by some systematic error - perhaps the wings of 
the Chandra psf at these energies may have been underestimated, or the novel method employed 
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by Woo et al to reconstruct the radial profile may be prone to systematic errors which are not fully 
understood. 

5.5. Discussion 

We have reviewed a number of measurements of dust scattering halos, and compared the 
predictions of the WD01 grain model to the values of T sca /Ay estimated from these observations. 
The results, shown in Fig. 11, are somewhat equivocal. Two studies (Woo et al. 1994; Predehl & 
Schmitt 1995) find values of T sca /Ay much lower than expected for the WD01 dust model, but we 
give arguments why these observations might have underestimated r sca . 

Draine & Tan (2003) have quantitatively modelled the X-ray halo around Nova Cygni 1992 
(typical photon energy ~0.48 keV) using the WD01 grain model. The halo intensity observed by 
ROSAT can be reproduced using a dust column density which agrees with the reddening inferred 
from the observed Ha/H/3 intensity ratio. 

Chandra observations of GX 13+1 by Smith et al. (2002), after correcting for missed halo 
counts using the constant dust gradient model with (3 = 1, give T sc& /Ay within a factor ~1.5 of the 
WD01 model at 2 keV, and in agreement at 3.5-4 keV. The corrections are somewhat sensitive to 
the (uncertain) dust density distribution for this case, since the halo intensity was not measured at 

e h < so". 

The recent Chandra measurement of the halo around Cyg X-l (Yao et al. 2003) implies values 
of T sca /Ay which are a factor ~1.5-2 smaller than expected for the WD01 model at E < 2.5 keV. 
At 3-4 keV the Yao et al. results are in excellent agreement with the WD01 model in terms of both 
TscaMv and the half-light radius of the halo. For 5-7 keV the values of T sca /Ay found by Yao et 
al. exceed the WD01 model by factors of 2-3, although the estimated errors are also large. 

Taken together, the Chandra observations of GX 13+1 (Smith et al. 2002) and of Cyg X-l 
(Yao et al. 2003) suggest that the WD01 model may have overestimated T sc& /Ay by a factor ~1.5 
between 1 and 2.5 keV. However, the Yao et al observations of Cyg X-l at 3-4 keV give both a halo 
concentration and r sc& /Ay in good agreement with the WD01 model, and the WD01 model is also 
consistent with the Nova Cygni observations at ~0.5 keV (Draine & Tan 2003). At this time we 
can conclude only that the WD01 model appears to give T sca /Ay in agreement with observations 
to within a factor ~1.5, but the existing observations do not permit a more precise statement. 

It is hoped that future observations by Chandra or XMM will be able to carry out high- 
signal-to-noise observations of X-ray scattering halos on sightlines where the dust distribution and 
reddening are well-determined. An optimal situation would be to use a source which is known to 
be distant compared to the dust doing the scattering, so that we can assume that x < 0.2 in eq. 
(12). An extragalactic X-ray point source (AGN or quasar) would be ideal for this purpose. 
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6. Summary 

The following are the principal results of this work: 

1. Dielectric functions for graphite and MgFeSi04 have been constructed which are continuous 
from submm to hard X-rays, obey the Kramers-Kronig relations, and satisfy the oscillator 
strength sum rule. 

2. Absorption, scattering, and extinction have been calculated for the WD01 grain model at 
X-ray energies. The calculated absorption edge structure appears to be consistent with re- 
cent spectroscopy by the Chandra X-Ray Observatory, and can be further tested by future 
observations. 

3. Differential scattering cross sections are presented for the Milky Way dust model at X- 
ray energies. These can be used for modelling X-ray scattering halos, and are available 
at http://www.astro.princeton.edu/~draine . 

4. The median scattering angle s ,5o(E) is given, as well as the scattering angles S) io and # S) go 
for 10% and 90% enclosed power. These can be used to assess the sensitivity of imaging 
observations for determination of the flux of halo photons. 

5. Simple analytic functions provide a good approximation to the differential scattering cross 
section for dust (eq. 10,11). 

6. We provide analytic approximations to g(0h), the fraction of the halo counts interior to Oh, for 
dust in a single sheet (eq. 17) and uniformly-distributed dust or dust with a density gradient 
(eq. 19). 

7. The total scattering cross section calculated for the WD01 grain model is compared with 
observations of X-ray halos by ASCA, ROSAT, and Chandra (see Fig. 11). The results are 
somewhat equivocal, and in some cases depend on corrections which are sensitive to the spatial 
distribution of the dust. ROSAT observations of Nova Cygni 1992 at ~0.5 keV (Draine k, 
Tan 2003), and Chandra observations at 3-4 keV of Cyg X-l (Yao et al. 2003) are in good 
agreement with the WD01 dust model, although the 2-3.4 keV observations of GX 13+1, and 
1-3 keV observations of Cyg X-l, give r sca /Ay about a factor ~1.5 below the WD01 model. 
At this time it is possible to conclude from the observations that r sca /^4y given by the WD01 
model is accurate to within a factor ~1.5, but a more precise statement is not possible. 

8. Analysis of the angular structure of X-ray halos around Galactic sources will generally be 
compromised by uncertainties concerning the location of the dust responsible for the scat- 
tering. The ideal observation is to observe an extragalactic point source, in which case the 
scattering dust is all at x « 0. It is hoped that such observations will be carried out by Chan- 
dra or XMM for bright extragalactic sources located behind sufficient Galactic dust, thereby 
providing a definitive test of this and other dust models. 
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